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ABSTRACT 

Laminar flow heat exchangers which cool oil in non-inter- 
connected parallel passages can experience non-uniform flows and 
a reduction in the effective heat exchanger coefficient in a range 
of Reynolds number which varies with tube length and diameter, 
tube wall temperature and fluid inlet temperature. The method of 
predicting the reduction in effective heat transfer coefficient 
and the range of Reynolds number over which these instabilities 
exist is presented for a particular oil, Mobil aviation oil 120. 
Included, also, is the prediction of the effect of radial viscosity 
variation on the constant property magnitudes of friction and heat 
transfer coefficient. 
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INTRODUCTION 

Heat exchangers that operate in the laminar flow regime with non- 
interconnected parallel passages can depart from predicted performance as 
a result on non-uniform flow distributions. Variations in the flow dis- 
tribution among the tubes can be caused by the effect of the fluid vis- 
cosity/temperature relation, natural convection and non-uniform passage 
sizes (from design and manufacturing tolerances) . If the exchanger is 
oriented vertically with the colder side up, the effect of natural con- 
vection can be eliminated. Non-uniform passage size can be minimized 
through quality control. Therefore, the effect of the viscosity/ tempera- 
ture relation on non-uniform flow will be analyzed in this paper. 

The viscosity/temperature relation for the Mobil aviation oil 

120 used in this analysis is presented in Fig. 1. Viscosity is a strong 

function of temperature, decreasing as temperature increases. In this 

analysis it is assumed that other fluid properties such as density and 

3 3 

thermal conductivity are uniform at p = 55.0 Ibm/ft (881.0 kg/m ), 
k = 0.08 BTU/hr-ff F (0.138 W/m*C) and Cp = 0.5 BTU/lbm-F (2.09 x 10^ 
J/kg*C ). This is a reasonable assumption for oils. Also, it is assumed 
that pressure drop entrance and exit effects can be neglected. The 
analysis was performed for flow in tubes. Uniform wall temperature was 
ussed assuming the heat transfer coefficient was very much greater on 
the outside than on the inside of the tube. 

When cooling liquids or heating gases, non-uniform flow can 
exist in the flow passages. For the same overall pressure drop, two 
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discrete flow rates can exist in the same heat exchanger. 

This problem has been discussed by Mueller [1], [2] and [3] 
and Rohsenow [4]. 

EFFECT OF RADIAL VISCOSITY VARIATION 

The fully developed fluid and temperature distributions in a 
tube are non-parabolic due to the viscosity/ temperature relation. This 
causes a lower velocity of the colder fluid near the tube wall and a 
higher velocity of the hotter fluid near the centerline compared with 
the constant property case (ie . uniform viscosity) with the same flow 
rate and bulk temperature. 

For steady state, fully developed flow and viscosity as a 
function of temperature, the momentum equation reduces to: 



dp _ 1_ ^ 
dx r 9r 



f 



( 1 ) 



Neglecting the effects of axial conduction and assuming constant heat 
flux at the tube wall, the temperature distribution for fully developed 
conditions is : 

An expression is needed to account for the change in friction 
factor and Nusselt number for the variable properties case (ie. variable 
viscosity) compared with the constant properties case (ie. uniform 
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viscosity) • Assume that the effect can be expressed as the ratio of the 
wall to bulk viscosities raised to an exponent as follows: 




f and Nu are based on constant properties at the bulk 

Cp Cp 

temperature . 

To find exponents m and n, Eqs. (1) and (2) were solved simul- 
taneously using an iterative procedure. Ah average fluid velocity, walL 
temperature and parabolic temperature distribution were assumed initially. 
The fluid viscosity versus temperature relation of Fig. 1 was closely 
approximated by a piecewise curve fit using a sixth degree polynomial. 
Based on the assumed temperature distribution, the viscosity |i(r) was 
determined. Equation (1) was solved for the fluid velocity based on the 
viscosity ]i(r). The dp/dx term was adjusted to obtain the desired 
average fluid velocity. Next, the temperature distribution from Eq. (2) 
was determined using the calculated fluid velocity. Then, the friction 
factor and Nusselt number were calculated. The viscosity |i(r) was de- 
termined for the revised temperature distribution and the procedure was 
repeated until the calculated friction factors and Nusselt numbers agreed 
within less than 1%. 

The results of this analysis showed that the exponents m and 



m 
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n varied as a function of wall and bulk temperature as shown in 
Figs. 2 and 3. The exponents were invariant to changes in the average 
fluid velocity. 

SINGLE TUBE PRESSURE DROP 

Having quantified the radial effects of the viscosity/tempera- 
ture relation on friction factor and Nusselt number using the ratio of 
the wall to bulk viscosities and the exponents m and n, an analysis 
is performed on a heat exchanger tube. An energy balance at the tube 
wall for uniform wall temperature yields: 



w c dT = hTT D(T - T )dx (5) 




Preliminary calculations showed that the fluid flow was fully 
developed a short distance from the tube entrance and entrance effects 
could be neglected. However, entrance effects on the temperature dis- 
tribution had to be taken into account because of the slowly developing 
temperature distribution. The following equation was suggested by 
Kays [5] to account for variation in the Nusselt number in the entrance 
region of tubes for uniform wall temperature; 



Nu = 3.66 + 



0.0668 TTx k 



1 + 



0.04 I ^ 

I TTx 



Wr 



2/3 



(7) 



10 



Therefore, the equation to determine the temperature drop along the heat 
exchanger is given by combining Eqs . (4), (6) and (7): 



AT 

T-T 

o 




wc 



0.0668 



TTx k 



1 + 0.04[ 



wc 



irx 



2/3) 



TTk 



WC 

p 



Ax 



( 8 ) 



Equation (8) includes the radial effects and the developing temperature 
distribution effects. The temperature change along the tube was solved 
by iteration of values of bulk viscosity, local temperature and exponent 
n . The temperature change and bulk viscosity along the tube for various 
flow rates is shown in Figs, 4 and 5 for a 0,25” diameter tube and tube 
wall temperature of 40 F, 

It can be shown that the pressure drop for an increment dx 
along a tube is given by: 



= W1g2 
dx Re D 2p 

Equation (9) is multiplied by Eq. (3) to account for radial effects on 
the friction factor. 



dp _ 64 
dx Re^ 




1 G 



D 2p 



( 10 ) 



In integral form Eq. (10) becomes: 
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( 11 ) 



The curves of pressure drop versus flow rate for various lengths 

of tube were calculated (Fig. 6) for D = 0.25", T = 40 F and T. = 300 F. 

o X 

It is observed that AP goes through a maximum, then a minimum, followed 
by a continuous rise as the flow rate is increased. From Eq. (11), 

AP - D ; therefore AP vs. w curves for other tube diameters would 
look like those of Fig. 6 but be displaced vertically proportional 
to D 



MULTI-TUBE HEAT EXCHANGER PERFORMANCE 

These curves are for single tubes but may be used to predict 
the flow maldistribution and resulting effective Nusselt number for 
multi-tube heat exchangers. This exchanger performance is determined 
by how the total flow is supplied - constant flow pump (gear pump), 
constant head source or variable head-flow pump (centrifugal pump). 

a) Constant Flow Pump 

In the case of a heat exchanger with five tubes (1 = 4^^, 

D = 0.25"), non-uniform flow can exist as shown in Fig. 7. For an in- 
let flow rate of 25. 75 Ibm/hr (average flow rate of 5.15 Ibm/hr per 
tube) and inlet temperature of 300 F and wall temperature of 40 F, 
four of the tubes have a flow rate of 1.26 Ibm/hr with an exit tempe- 
rature of 42.7 F and one tube has a higher flow rate of 20.7 Ibm/hr 
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with an exit temperature of 220.0 F. The average bulk temperature at 
the exchanger outlet is 185.2 F. Figure 8 shows a plot of temperature 
along the exchanger for one tube at higher flow and four tubes with 
lower flow. The average exit temperature (T^^g) is the summation of 
the individual tube exit temperature weighted by the percentage of 
flow through each tube. This is the temperature that would be mea- 
sured at the heat exchanger exit after mixing the five streams. If 
it were not known that this maldistribution existed and it were as- 
sumed that the flow was uniform in the five tubes at one-fifth of 
the total flow, the effective Nusselt number would be tetermined 
using T^vg follows: 




( 12 ) 



For this case ~ 5.72 BTU/hr *f t^*F (32.49 W/m^*C) and ” 1*49. 

Also shown in Fig. 8 is the temperature along the exchanger and the exit 

temperature (T ) if the tubes, in fact, had uniform flow. For 

unirorm 

this case h .. = 10.37 BTU/hr-ft^-F (58.87 W/m^*C) and Nu .. = 2.70. 
err etr 

Since T > T , the Nusselt number (and heat exchanger per- 

avg uniform 

fo nuance) is lower when non-uniform flows exist. 

An enlarged view of the region of non-uniform flows is shown 
in Fig. 9. On Fig. 9 curve (c-d) at any AP level represents the 
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average flow rate with four tubes having a flow rate on curve (u-a) and 

one tube on curve (m-1) , the stable flow curves. Then 5w ^ = 4w H- w « , 

cd ua roio 

Curves (s-f) , (g-h) and (o-j) are similar curves for flow splits of 3-2, 

2-3 and 1-4 tubes on curves (u-a) and (m-1) . As the total flow rate is 

increased, uniform flow exists through the peak at point (a) and on to 

point (b). At point (b) the pressure drop decreases suddenly to point 

(c) where one tube has a higher flow rate while the other four tubes 

have lower flow rates. The total flow rate is: 5w- = 5w = 4w + w . 

b c u m 

As flow is further increased, flow increases in each of the tubes 
(along the pressure drop curve) with the average flow rate shown as 
line (c-d) . Point (d) is the point depicted in Figs. 7 and 8. As flow 
is increased slightly above the flow corresponding to point (d) , the 
pressure drop falls from points (d) to (e) and now two tubes have 
.higher flow rates and three tubes have lower flow rates. As the flow 
is increased, the flow in each tube increases with the average flow 
rate shown on the line (e-f ) . This process continues until uniform 
flow in all five tubes occurs at point (k) . In the case of infinite 
tubes, the average flow rate would follow the line (a-1) with flow in 
the tubes corresponding to the flow at either point (a) - (lower flow 
rate) - (1) - (higher flow rate). The mix of tubes with higher and 
lower flow rates is determined by the average fluid flow rate. 

For the case of decreasing flow rate, non-uniform flows also 
occur. As flow is decreased to point (m) the pressure drop decreases 
with all tubes having uniform flow. The pressure drop suddenly in- 
creases from points (m) to (n) at which one tube has lower flow and 
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four tubes have higher flow. As flow is decreased further, flow in each 
tube decreases with the average flow rate given by line (n-o). Addi- 
tional tubes "pop-over" to lower flow at points (o), (q), (s) and (c). 
All the tubes have the same flow rate at point (b). In the case of 
infinite tubes the average flow rate would follow line (m-u) with flow 
in the tubes corresponding to the flow at either point (u) - lower 
flow rate - or (m) - higher flow rate. 

Figure 10 shows the effects of non-uniform flow on the Nusselt 
number as a function of Reynolds number for the case of five tubes. 

The solid line is for increasing flow and the dashed line is for de- 
creasing flow. The Nusselt number is reduced significantly between 

points (b) and (k) for increasing flow and between points (m) and (b) 

% 

for decreasing flow. These are the regions of the curve where non- 
uniform flow exists. 

For the case of infinite tubes, the effect of non-uniform flow 
on exchanger performance is shown in Figs. 11 and 12. As the tube 
diameter is decreased, the region of non-uniform flow is extended 
to higher Reynolds numbers but the Nusselt number decreases less. 
Increasing the tube length has an effect on exchanger performance 
similar to that of decreasing the tube diameter. 

b) Constant Head Source 

For the constant head source, as the pressure drop is in- 
creased beyond point (a) of Fig. 9, the average flow rate will in- 
crease to point (1) regardless of the number of tubes. With con- 
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tinual increase in constant head, operation can not be sustained at 
flow rates between points (a) and (1). For the case of continual de- 
creasing constant head, all tubes remain operating along line (1) to 
(m) . Then the flow drops to point (u) in all tubes. 

c) Variable Head-Flow Pump 

For a centrifugal pump source with the pump characteristic 
shown as lines (a-b) and (c-d) in Fig. 13, non-uniform flows can also 
exist. The situation is much the same as for the case of a constant 
head source except that the average flow follows the pump characteristic 
instead of the constant head line. As pump speed is increased above 
the point (a) all tubes shift to operation at point (b). The ex- 
changer can not operate between points (a) and (b) at the higher 
speed nor between points (c) and (d) at the lower speed. 



DISCUSSION 

Variations in the tube wall and inlet temperature can affect 
the values of exponents ra and n (Figs. 2 and 3). Together with 
the y vs. T relation (Fig. 1), changes in exponents m and n 
can cause differences in the shape of the AP vs. w curve determined 
by Eq. (11). The AP vs. w curves for a single tube (L = 4^, 

D = 0.25**) are plotted in Fig. 14 for different combinations of wall 
and inlet temperature. The tube with T^ = 40®F and T^ = 300°F 
exhibits a region in which non-uniform flow can exist. If T^ is 
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changed from 300°F to 200®F, there is only a slight dip in the curve 
and the region of non-uniform flow is narrower. If is changed 

from 40®F to 100®F with = 300®F, no region of non-uniform flow 
exis ts . 

Results were similar for calculations done using another oil, 
Mobil automotive oil HD-30, 

CONCLUSIONS 

Laminar flow, parallel passage heat exchangers can experience 
maldistribution of flow and a reduction in the effective heat transfer 
coefficient in a range of Reynolds number which varies with length, 
tube diameter and the levels of the temperature - wall and inlet, 

For Mobil aviation oil 120 with the y vs, T relation of 
Fig, 1, the exponents m and n in the viscosity correction ratios, 
Eqs, (3) and (4), are shown to vary with both wall and bulk temperatures. 
Figs. 2 and 3. 
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NOMENCLATURE 



c 

P 



D 




G 

k 

L 



m 

n 

n 

o 

^"eff 

Nu 

Nup 

P 

AP 



T 

T 

avg 

T 



specific heat at constant pressure 
tube diameter 

friction factor (variable viscosity) 

friction factor (uniform viscosity evaluated at the bulk 

temperature) 

mass flow rate 

effective heat transfer coefficient 
thermal conductivity 
tube length 

exponent for friction factor correction 
exponent for Nusselt number correction 
number of tubes ^ 

effective Nusselt number (variable viscosity) 

Nusselt number (variable viscosity) 

Nusselt number (uniform viscosity at the bulk temperature) 

pressure along the tube 

tube pressure drop 

tube radius, radial direction 

Reynolds number evaluated at \i^ 

temperature as function of r or x 

average exit temperature in non-uniform flow condition 
bulk temperature 
wall temperature 
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inlet temperature 

"^uniform average exit temperature for uniform flow condition 
u fluid velocity in x - direction 

W total flow rate 

w individual tube flow rate 

X direction along tube 

Greek Symbols 



P 

y 

^^0 



% 



fluid density 

viscosity as function of T 
viscosity evaluated at 
viscosity evaluated at T, 

D 
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FIGURE I VISCOSITY VS TEMPERATURE RELATION 
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FIGURE 4 TEMPERATURE VARIATION ALONG EXCHANGER 



BULK VISCOSITY, Ibm/hr.ft 
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FIGURE 5 VISCOSITY VARIATION ALONG EXCHANGER 



PRESSURE DROP, Ibf /ft 
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TEMPERATURE, 



26 




FIGURE 7 FIVE TUBE HEAT EXCHANGER 




FIGURE 8 TEMPERATURE DISTRIBUTION FOR FIVE TUBE EXCHANGER 
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FIGURE 9 FLOW MALDISTRIBUTION IN FIVE TUBE EXCHANGER 



EFFECTIVE NUSSELT NUMBER 
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FIGURE 10 FIVE TUBE EXCHANGER PERFORMANCE 
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FIGURE II EFFECT OF DIAMETER ON PERFORMANCE 
FOR AN INFINITE NUMBER OF TUBES 



EFFECTIVE NUSSELT NUMBER 
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FIGURE 12 EFFECT OF LENGTH ON PERFORMANCE FOR AN INFINITE 
NUMBER OF TUBES 



PRESSURE DROP. Ib|/tt 
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FIGURE 13 AP VERSUS W FOR VARIABLE HEAD-FLOW 
PUMP FOR INFINITE NUMBER OF TUBES 



PRESSURE DROP, Ibf/ft 
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FIGURE 14 EFFECT OF TEMPERATURE LEVELS ON PRESSURE 
DROP CURVES 
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